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Abstract 

Conformo static charged dust distributions are investigated in the 
framework of general relativity. Einstein's equations reduce to a non- 
linear version of Poisson's equation and Maxwell's equations imply the 
equality of the charge and mass densities. An interior solution to the 
extreme Reissner-Nordstrom metric is given. Dust distributions con- 
centrated on regular surfaces are discussed and a complete solution is 
given for a spherical thin shell. 







1 Introduction 



Let M be a four dimensional spacetime with the metric 

g^v = f^rii^v-u^Uv (1) 

where r\ liliV = diag (0, 1, 1, 1) and = v7^°- Here Latin letters represent 
the space indices and 8y is the three dimensional Kronecker delta. In this 
work we shall use the same convention as in [I]. The only difference is that 
we use Greek letters for four dimensional indices. Here M is static. The 
inverse metric is given by 

sr = fi£ v -vru" (2) 
where 7? 2 ^ = diag (0, 1, 1, 1) , u M = g fJ,u u u = — h= 5q. Here is a time-like 
four vector, u^u^ = —1. 

The Maxwell antisymmetric tensor and the corresponding energy momentum 
tensor are respectively given by 

F llv = V v A ll -V li A v (3) 
M, v = ±-{F, a F:-- A F 2 g, v ) (4) 

where F 2 = F^ v F fJiU . The current vector j M is defined as 

V u F» v = Aixf (5) 
The Einstein field equations for a charged dust distribution are given by 

G flu = 8ttT^ u = + (8tt p)u fl u v (6) 

where p is the energy density of the charged dust distribution and the four 
velocity of the dust is the same vector appearing in the metric tensor. 
Very recently we invetstigated the above filed equations with metric given 
in flU). We find that j M = p e M M , where p e is the charge density of the dust 
distribution. Let A be a real function depending on the space coordinates. 
It turns out that A4 = and 
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where = ±1. Then the field equations reduce simply to the following 
equations. 



V 2 A + 47rpA 3 = (8) 
p e = kp (9) 

where V 2 denotes the three dimensional Laplace operator in Cartesian flat 
coordinates. These equations represent the Einstein and Maxwell equations 
respectively. In particular the first equation ([8]) is a generalization of the Pois- 
son's potential equation in Newtonian gravity. When p vanishes , the space- 
time metric describes the Majumdar-Papapetrou space-times ' @ ' @] ' [§]• 
For the case p ^ , the reduced form of the field equations (H) were given 
quite recently (see also [|J). 



2 Charged dust clouds 

In the Newtonian approximation X = 1 + V , Eq.(||) reduces to the Poisson 
equation, V 2 V + 4tt p = 0. Hence for any physical mass density p of the dust 
distribution we solve the equation (§) to find the function A. This determines 
the space-time metric completely. As an example for a constant mass density 
P = Po > we find that 

\ = -^ = cn(l l x i ) (10) 

Here U is a constant three vector ,a 2 = 1$ and cn is one of the Jacobi elliptic 
function with modulus square equals |. This is a model universe which is 
filled by a (extreme) charged dust with a constant mass density. 



3 Interior solutions 

In an asymptotically flat space-time , the function A asymptotically obeys 
the boundary condition A —>■ Aq (a constant). In this case we can establish 
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the equality of mass and charge e = ±mo, where mo = / p \/—g d 3 x. For 
physical considerations our extended MP space-times may be divided into 
inner and outer regions. The interior and outer regions are defined as the 
regions where Pi > and p = respectively. Here % = 1,2,...,N, where N 
represents the number of regions. The gravitational fields of the outer regions 
are described by any solution of the Laplace equation V 2 A = , for instance 
by the MP metrics. As an example the extreme Reissner-Nordstrom (RN) 

metric (for r > R ), A = Ao H — - may be matched to a metric with 

r 

sin(6r) , . 

A = a — - — -, r < Ro 11 

r 

describing the gravitational field of an inner region filled by a spherically 
symmetric charged dust distribution with a mass density 

P = P(0)[-^] 2 (12) 
sin (or) 

Here p(0) = t-^t, r 2 = Xi x\ a and b are constants to be determined in terms 
of the radius Rq of the boundary and total mass mo (or in terms of p(0)). 
The boundary condition , when reduced on the function A on the surface 
r = R it must satisfy both \ out = X in and \' out = X' in . Here prime denotes 
differentiation with respect r. They lead to 



m = ^, (13) 

Ao = abcos(bRo), (14) 

Ai = a (sin(6_R ) — bR cos(6_R )) (15) 

When the coordinates transformed to the Schwarzschild coordinates .i.e., 
Ao r + Ai —>■ r then the line element becomes 



ds z=-^ { l-±fdt* + -^ + r>dtf. (10) 

v r ) 



hence Ai is the mass in the Newtonian approximation then Ai = rriQ. 

In this way one may eliminate the singularities of the outer solutions by 

matching them to an inner solution with a physical mass density. 
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For the mass density p = in general , we may have the complete 
solution. Here b is a nonzero constant which is related to mo by the relation 
b 2 Rq = 3m and we find that 

A = £a<,mjz(&0*WM) (17) 

l,m 

where ji(br) are the spherical Bessel functions which are given by 

and Y^ m are the spherical harmonics. The constants a^ m are determined 
when this solution is matched to an outer solution with V 2 A = 0. The 
interior solution given above for the extreme RN metric with density (|Tj) 
corresponds to I = 0. 



4 Point particle solutions 

Newtonian gravitation is governed by the Poisson type of linear equation. 
Gravitational fields of spherical objects in the exterior regions may also be 
identified as the gravitational fields of masses located at a discrete points 
(point particles located at centers of the spheres) in space (i? 3 ). The solution 
of the Poisson equation with N point singularities may be given by 

N m- 

A = l + E^ (19) 

i=l '* 

n = [(x - Xi) 2 + (y - Vl ) 2 + {z- Zl ) 2 } 3 (20) 

where N point particles with masses m 8 are located at the points (x l ,y l , z % ) 
with i = 1,2, ,,N. The same solution given above may also describe the 
exterior solution of the iV spherical objects (with nonempty interiors) with 
total masses rrii and radii R{. The interior gravitational fields of such spher- 
ical objects can be determined when the mass densities pi are given. The 
essential point here is that the limit Ri — > is allowed. This means that 
the dust distribution is replaced by a distribution concentrated at the points 
(x l , y % .,z l ) . Namely in this limit mass densities behave as the Dirac delta 
functions, 
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N 

p ^J2mi5(x - Xi)8(y - Ui)S(z - Zi). 

i=l 

This is consistent with the Poisson equation V 2 A + Air p = , because — in 
the solution (|19| ) is the Green's function, i.e., 

V 2 — = -47r 5(x - Xi) 5(y- y») 5(z - z t ). 

Such a limit , i.e., Ri — > is not consistent in our case , V 2 A + 47r p A 3 = 0. 
The potential equation is nonlinear and in particular in this limit the product 
of p and A 3 does not make sense. Hence we remark that the Majumdar- 
Papapetrou metrics should represent the gravitational field N objects with 
nonempty interiors (not point-like objects). 

5 Thin shell solutions 

In the previous section we concluded that the dust distribution can not be 
concentrated to a point. We observe that , the potential equation (§) does not 
also admit dust distributions on one dimensional (string like distributions) 
structures. This is compatible with the results of Geroch-Traschen ||. On 
the other hand , the mass distribution p can be defined on surfaces. 

Let S be a regular surface in space (R 3 ) defined by S = [(x, y, z) G 
R 3 ; F(x, y, z) = 0] , where F is a differentiable function in R 3 . When the 
dust distribution is concentrated on S the mass density may be represented 
by the Dirac delta function 

p(x, y, z) = p (x, y, z) 5(F(x, y, z)) (21) 

where p (x,y,z) is a function of (x,y,z) which is defined on S. The func- 
tion A satisfying the potential equation (§) compatible with such shell like 
distributions may given as 

A(x, y, z) = A (x, y, z) - \ x (x, y, z) 9(F) (22) 

where Ao and Ai are differentiable functions of (x, y, z) and 0(F) is the Heav- 
iside step function. With these assumptions we obtain 
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1 VAx-VF , 

Po(x, y,z) = — j—^ | s (23) 

V 2 A = V 2 Ai = (24) 
and in addition Ai|s = 0. We have some examples: 

1. S is the plane z = 0. We have p(x,y,z) = p (x,y) S(z). Then it follows 
that 



Hx,y,z) = X (x,y,z) - X 2 (x,y) z9(z) (25) 

Po(w) = iix^F (26) 

V 2 A = V 2 A 2 = (27) 



2. S is the cylinder F = r — a = 0. We have p(r, 8, z) = p (9, z) 5(r — a). 
Then it follows that 



X(r,9,z) = X (r,e,z) - X 2 (9,z) \n(r/a)6(p - a) (28) 

Po(r,9,z) = t^tt^ (29) 
Ana (Xq\ s ) 6 

V 2 A = V 2 A 2 = (30) 

Here we remark that the limit a — > does not exist. This means that the 
mass distribution on the whole z- axes is not allowed. 

3. S is the sphere F = r — a = 0. We have p(r, 9, 0) = Po(#> 0) <K r — a )- 
Then it follows that 



X(r,9,<f>) = Xo{r,e,4>)-X 2 {e,4>){---)e{r-a) (31) 

a r 

= 4^ (W <32) 

V 2 A = V 2 A 2 = (33) 
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We note that the total mass is infinite on non-compact surfaces. 

For compact case we shall consider the sphere in more detail. In this case 
we may have A = fi A 2 + if) such that V 2 ip = and if) (a, 9, 0) = 0. We shall 
assume if> = everywhere, then 

Po = 4^ (34) 

Hence the total mass mo on S is given by mo = / sj—g pd 3 x = -. Let 
yout anc j yin ojgj^ote solutions of (H) corresponding to the exterior and inner 
regions respectively. They are given by 



X out (r,9,(P) = A(r>a,#,0) = l-^ + ^ (35) 

a r 

\ in (r, 9, 0) = A(r < a, 9, 0) = 1, (36) 

where we let A2 = mo- Here we remark that the point particle limit a — > 
does not exist. The solution given above represents the extreme Reissner- 
Nordstrom solution 



ds 2 = - (1 - ^ + ^)~ 2 dt 2 + (1 - ^ + ^) 2 (dr 2 + r 2 d n 2 ) (37) 
a r a r 

By letting r = R ~™° where f3 = 1 — ^ and a 7^ m is assumed. We obtain 
the extreme RN in its usual form 

^ = - P 2 (1 - ^) 2 ^ 2 + ^^y 2 + R 2 dQ 2 (38) 

Hence we obtain a solution where the exterior solution is the extreme Reissner 
- Nordstrom metric, but inside the sphere with radius a, the spacetime is flat. 
Thus extreme RN solution is matched to a spherical shell (R = a) of dust 
distrubition. 

The case a = mo represents the Levi-Civita -Bertotti-Robinson spacetime 
outside the dust shell and the flat spacetime inside. 

ds 2 = -—, dt 2 + H$dr 2 + m 2 d n 2 (39) 
mg r 2 
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By letting r = we obtain the usual conformally flat LCBR metric 

ds 2 = ^§ [-dt 2 + rfi? 2 + R 2 dtt 2 } (40) 

In the new coordinates the surface is again R = mo. In these dust shell 
solutions the function A is continuous on the surface r = a , but its normal 
derivative to S is discontinuous , 

Kut - Kn = -47r(T = -m /a 2 
as expected, where a = mass per unit area = m /47ra 2 . For thin shells in 



general relativity see the recent work of Mansouri and Khorrami |10| and 
also Mansouri's contribution in this proceeding. 



6 Conclusion 

We have solved the Einstein field equations in a conformo-static space-time 
for a charged dust distribution. We reduced the whole Einstein field equa- 
tions to a nonlinear Poisson type of potential equation (||). Physically rea- 
sonable solutions of this equation give an interior solution to an exterior MP 
metrics. We have given some explicit exact solutions corresponding to some 
mass densities. In particular we have given an interior solution of the extreme 
RN metric. 

We showed that the limiting cases of mass distributions on discrete points 
and also on lines in R 3 are not possible. We have examined some possible 
mass distributions on regular surfaces. We have found solutions correspond- 
ing to shell like dust distributions. In particular for the spherical dust shell 
we presented an exact solution given in eqns. fl35[ - |39| ) representing a spherical 
shell cavity immersed in the extreme RN spacetime. 
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